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The elastic constant tensors for the hep phases of three transition metals (Co, Re, and Fe) are 
computed as functions of pressure using the Linearized Augmented Plane Wave method with both 
the local density and generalized gradient approximations. Spin-polarized states are found to be 
. stable for Co (ferromagnetic) and Fe (antiferromagnetic at low pressure). The elastic constants 

' of Co and Re are compared to experimental measurements near ambient conditions and excellent 

agreement is found. Recent measurements of the lattice strain in high pressure experiments when 
interpreted in terms of elastic constants for Re and Fe are inconsistent with the calculated moduli. 

o 
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I. INTRODUCTION 



The effect of pressure on the propagation of elastic waves in materials is essential for understanding interatomic 
interactions, mechanical stability of solids, phase transition mechanisms, material strength, and the internal structure 
of the Earth and other planets. However, little is known of the elasticity of solids at high pressure. The experimental 
study of the elasticity of materials under high pressure is challenging, as traditional methods have been applied only 
to moderate pressures. Ultrasonic measurements are generally limited to a few GPaLU, while Brillouin spectroscopy 
has been applied up to 25 GPaB. 

We investigate the elasticity of three hexagonal transition metals at high pressure: iron, rhenium, and cobalt. 
High pressure properties of iron are of considerable geophysical interest as the Earth's solid inner core is composed 
priraadly of this element. The elasticity of hep iron is important for understanding the elastic anisotropy of the inner 
■ coretTEl, and its super-rotationD. Rhenium is the strongest metal known at high pressured and is widely used as a 
gasket material in diamond anvil cell experiments. We have chosen cobalt for this study because of its proximity to 
■^j- . iron in the periodic table and as an example of a ferromagnetic hep metal. 

All three of these metals have been studied experimentally under high pressure and thfiir equations-of-state are well 
0^ , known. Iron transforms from the bec phase.-at ambient conditions to hep near 13 GPaQ; the equation-of-state of the 
hep phase has been measured up to 300 GPaO. Recent advances in diamond anvil cell techniques have made it possible 
to evaluate the lattice strain in a polycrystal subjected to a non-hydrostatic stress field which can be associated with 
elastic constants. The elasticity of iron has been inferred bv-this method at high pressure (up to 210 GPa)li!JO. The 
equation-of-state of cobalt has been measured up to 80 GPall3 and its elastic constants were obtained at pro pressure 
. using traditional ultrasonic methodstil In the case of rhenium the equation-of-state is know-u to 215 GPalij, its elastic 
constants and their pressure derivatives have been ultrasonically measured at low pressuret-3. The same experimental 
method for evaluating lattice strains as in hep iron has been applied to rhenium in the pressure range 18-37 GPali3. 
O ■ Iron has been studied widely with first-principles theoretical approaches because of its geophysical importance and 
the well known failure of the local .density approximation (LDA) to the exchange-correlation potential to predict 
the ferromagnetic bec ground statrP- 7 ! This failure was a major impetus in the development of the generalized 
gradient aoutpximation (GGA)E5~l3. The equation-of-state of hep iron under LDA and GGA is well known to high 
pressuresclTcfLand its elastic constants have been calculated bv_the full-potential linearized muffin-tin orbital method 
(FP-LMTO)Ej, and a total energy tight-binding (TB) methoaHO For hep cobalt calculations have been performed 
with the LMTO method in the atomic sphere approximation for LDAc3 and the linearized combination of atomic 
orbital method (LCAO) for GGAED. There is no previous theoretical woxk on the elastic constants of hep cobalt. For 
rhenium only one study has focused on the hep phase at high pressures; using FP-LMTO with LDA the equation- 
of-state and the elastic constants at zero pressure have been calculated. 

We organize the paper as follows. Section || elaborates the computational details of our first principles calculations 
and our approach to calculating the elastic constants, the elastic wave velocities, and the acoustic anisotropy. It 
is followed by a section presenting our results on the magnetic state of the materials studied, their c/a ratios, the 
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equation-of-state, and the elastic constants as functions of pressure. We compare our results in terms of the elastic 
wave velocities to high pressure experiments and the Earth's inner core. In section IV we analyze the elastic anisotropy 
resulting from our calculations, recent experimental and theoretical results, and the predictions from a central nearest 
neighbor force model. Finally, we present our conclusions in section |v|. 



II. METHOD 



LAPW Total Energy Calculations 

We investigate the energetics of hep iron, cobalt, and rhenium using the full-potential linearized-augmented plane- 
wave method (LAPW)E3 with-both LDA and GGA approximations to the exchange-correlation potential. For LDA the 
form of Hedin and LundquistEJ and von Barth and HedinEH are used for non-magnetic and spin-polarized calculations, 
respectively. For GGA we adopt the efficient formulation of Perdew, Burke, and ErnzhofcrE3. 

Core states are treated self-consistently using the full Dirac equation for the spherical part of the potential, while 
valence states are treated in a semirelativistic approximation neglecting spin-orbit coupling. We investigate ferromag- 
netic alignment in spin-polarized calculations for all metals and antiferromagnetism for iron. For consistency of the 
results all parameters in the calculations except for spin-polarization are kept fixed. 

For the 3d metals 3s, 3p, 3d, 4s, and Ap states are treated as valence electrons for all volumes. For rhenium we 
treat all electrons up to 4/ as core, 5d and 6s as valence states. For rhenium we also have tested other configurations, 
such as including the 4/ as valence states, which did not change our results significantly. The muffin-tin radii Rmt 
are 2.0 Bohr for the 3d metals, and 2.3 Bohr for rhenium. As spin-orbit coupling of the valence electrons is important 
for the band structure and other properties of heavy elements, we consider the influence of the spin-orbit term on the 
equation-of-state for Re by including it in a variational stepcE 

We carry out total energy calculations over a wide range of volumes for all three metals (0.7 Vq - 1.2 Vq, with 
Vo the zero pressure volume). At each volume we determine the equilibrium ratio of the lattice constants c/a by 
performing calculations for several different values of this ratio. The equatioiuof-state is obtained by describing the 
energy-volume curve with a third order expansion in the Eulerian finite strainE3. 

We have performed convergence tests with respect to Brillouin zone sampling and the size of the basis set, 
RMTK max , where K max is the largest-reciprocal space wave- vector in the basis set. Converged results are achieved 
with a 12x12x12 special k-point mesht^l, yielding 114 k-points in the irreducible wedge of the Brillouin zone for the 
hep lattice, and up to 468 k-points for the monoclinic lattice used in elastic constants calculations. The number of 
k-points in the. full Brillouin zone is well above the convergence criterion for elastic constant calculations established 
by Fast et alB3. The size of the basis set is given by RuT^max = 9.0, yielding 158 to 311 basis functions, depending 
on volume. Careful convergence tests show that with these parameters relative energies are converged to better than 
0.1 mRy/atom, magnetic moments to better than 0.05 /x^/atom, and c/a to within 0.025. 



Elastic Constants 



We calculate the elastic constants as the second derivatives of the internal energy with respect to the strain tensor 
(e). We choose the applied strains to be isochoric (volume-conserving) which has several important consequences: 
First, we assure the identity of our calculated elastic constants with the stress-strain coefficients, which-are appropriate 
for the calculation of elastic wave velocities; this identity is non-trivial for finite applied pressureoa. Second, the 
total energy depends on the volume much more strongly than on strain; by choosing volume conserving strains we 
obviate the separation of these two contributions to the total energy. Third, the change in the basis set associated 
with the applied strain is minimized, thereby minimizing computational uncertainties. 

We obtain the elastic constants at the equilibrium relaxed structure at any volume V by straining the lattice, 
relaxing the symmetry allowed internal degrees of freedom, and evaluating the total energy changes due to the strain 
as a function of its magnitude S. 

The bulk modulus K is calculated by differentiating the equation-of-state. For hexagonal crystals K is the combi- 
nation of elastic constants 

K = [C 33 (Cn + Ci 2 ) - 2C? 3 ] fC s , (1) 

with 

C s = C u + C12 + 2C 33 - 4C 13 . (2) 



2 



The volume dependence of the optimized c/a is related to the difference in the linear compressibilities along the a- 
and c-axes (k a and k c ). The dimensionless quantity R describes this as 

In terms of the elastic constants, 

R=(C 33 -C n -C 12 + C 13 )/C s . (4) 

We calculate Cs by varying the c/a ratio at a given volume, according to the isochoric strain 

(50 \ 

e{5) =0<50. (5) 
V (l + 5y 2 -lj 

The corresponding energy change is 

E{5) = E(0) + C S V5 2 + 0(5 3 ). (6) 

In the expressions for Cs, K, and R, C\\ and C\i occur only as a sum. To separate these constants we determine 
their difference, Cxi — C12 = 2Cqq by applying an orthorhombic strain, space group Cmcm. For the strained lattice we 
use the two atom primitive unit cell, with the atoms in the Wyckoff position 4c, coordinates (y,—y,l/4). The strain is 

(50 \ 

s(5) = -5 (7) 

\ 5 2 /(l-5 2 ) J 

leading to a change in total energy: 

E(5) = E(0) + 2C 66 V5 2 + O(S i ). (8) 

In the unstrained lattice the atomic coordinate is y — 1/3, but varies under strain!!. We relax our calculations with 
respect to this internal degree of freedom. 

To determine C44 we use a monoclinic strain, space group C2/m. The atomic positions in the two atom primitive 
unit cell are (1/6,5/6, 1/4). The strain applied 

/ 5 \ 

e(5) = 5 2 / (1 - 5 2 ) (9) 
\ 5 / 

results in an energy change 

E(5) = E(0)+2C U V5 2 +0(5 4 ). (10) 

The equilibrium positions of the atoms are unaffected by this strain and do not need to be redeterminedS 

While for C§q and C44 the leading error term is of the order S 4 , for Cs it is of third order in 5. It is therefore crucial 
to include positive and negative strains in the calculation for Cs- The strain amplitudes applied are typically nine 
values of 5 covering ± 4% for Cs; for Cqq and C44, seven values of 5 ranging to 6% are applied. The elastic constants 
are then given by the quadratic coefficient of polynomial fits to the total energy results; the order of the polynomial 
fit is determined by a method outlined by Mehltll. 

From— the full elastic constant tensor we can determine the shear modulus \i according to the Voigt-Reuss-Hill 
schemeES and hence the isotropically averaged aggregate velocities for compressional (v p ) and shear waves (v s ) 

, V s = \ZJjJp, (11) 

with p the density. 

More generally, the acoustic velocities are related to the elastic constants by the Christoffcl equation 

{CijkinjHk - M5u) Ui = 0, (12) 
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where Cijki is the fourth rank tensor description of elastic constants, n is the propagation direction, u the polarization 
vector, M = pv 2 is the modulus of propagation and v the velocity. 
The acoustic anisotropy can be described as 

_ Mi[n x ] 

1 = MpO]' (13) 

where n x is the extremal propagation direction other than [100] and i is the index for the three types of elastic waves 
(one longitudinal and the two polarizations of the shear wave). Solving the Cristoffel equation for the hexagonal 
lattice one can calculate the anisotropy of the compressional (P) wave as 

Ap = ^- (14) 
tu 

For the shear waves the wave polarized perpendicular to the basal plane (SI) and the one polarized in the basal plane 
(52) have the anisotropics 

. C\\ + C33 — 2C13 a C44 

Asi = — , A S2 = 7^-. (15) 

4O44 O66 

While for S2- and P-waves the extremum occurs along the c-axis, for SI it is at an angle of 45° from the c-axis in the 
a-c-plane. We note that an additional extremum may occur for the compressional wave propagation at intermediate 
directions depending on the values of the elastic constants. 

III. RESULTS 
Magnetism 

We find a stable ferromagnetic state only in cobalt. It is stabilized over a wide volume range with the magnitude 
of the moment decreasing with pressure in agreement with previous theoretical results on the pressure dependence 
of magnetic momentstJ in other transition metals. Only at the smallest volume considered (50 Bohr 3 , 180 GPa) is 
the moment vanishingly small (Fig. [I]). LDA and GGA yield consistejut results and predict a zero pressure magnetic 
moment of 1.55 fj,B, in excellent agreement with experiment (1.58 /xblj). 

In the case of hep iron, we also investigate two antiferromagnetic states. The first consists of atomic layers of 
opposing spin perpendicular to the c-axis (afml). The other arranges the planes of opposite spins normal to the [100] 
direction in the hep lattice; this can be described by the orthorhombic representation of the hep unit cell (space group 
Pmma) with spin up in the (1/4,0,1/3) and spin down in the (1/4,1/2,5/6) position (afmll). We find that both 
structures are more stable than the non spin-polarized state and that afmll is energetically favored over afml. For 
both antiferromagnetic states the moment is strongly pressure dependent ■— For afml it vanishes at volumes larger than 
Vq (Fig. [I]), in excellent agreement with results of Asada and TerakuraE3. The other structure, afmll, possesses a 
magnetic moment well into the stable pressure regime of hep iron, up to ~ 40 GPa. (Fig. [j]). Because of frustration 
on the triapgular lattice, it is possible that more complex spin arrangements such as incommensurate spin waves as 
for fee ironEil or a spin glass are still more energetically favorable .ihan afmll. 

Diamond anvil cell in situ Mossbauer measurements of hep iroro have shown no evidence of magnetism in the hep 
phase. The low antiferromagnetic moment we calculate in the stable hep regime and the significant hysteresis of the 
bec-hep transitioned might explain that no magnetism in hep iron has been detected in the high pressure Mossbauer 
experiment. In this context it may be relevant that indirect evidence for magnetism exists at-low pressure. Epitaxally 
grown iron-ruthenium, supcrlattices have shown magnetism occurring in hep iron multilayersc3. Its character, however, 
is still controversialtElcfl. 

c/a Ratios 

For all materials studied the c/a ratio agrees with experimental data to within 2% and is essentially independent 
of the exchange correlation potential (GGA or LDA). Equilibrium c/a ratios_for iron range from 1.58 at zero pressure 
to 1.595 at 320 GPa. This is consistent with experimental measurementstm in the range of 15 to 300 GPa, which 
have shown considerable scatter. For cobalt the zero pressure c/a ratio is calculated as 1.615, increasing to 1.62 at a 
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pressure of almost 200 GPa. The zero pressure c/a is slightly lower than the experimental ralue of 1.623E3. Diamond 
anvil cell experiments have found a higher value of c/a, as much as the ideal value (1.633)113, this discrepancy might 
be due to the coexistence of hep and metastable fee cobalt in the polycrystalline sampldla. The c/a ratio for rhenium 
(1.615)|-does not change over the whole pressure range studied - and is in good agreement with experimental results 

(i.6i3)ta 



Equation-of-State 

For the equation-of-state of rhenium, LDA shows better agreement with experimental data than does GGA (Fig. 
^, Table ||). GGA overestimates the zero pressure volume and softens the b ulk m odulus, supporting a general pattern 
seen in prior density functional calculations using GGA for other 5d metalsErEJ. Including spin-orbit coupling in the 
calculation has little effect on the equation-of-state parameters, resulting in less than 1% change in the zero pressure 
volume and 2% in the bulk modulus. For cobalt, as for other 3d metals GGA is superior to LDA and reproduces the 
experimental equation-of-state to within 2% in volume and 10% in bulk modulus (Fig. ||, Table |). 

The discrepancy in the equation-of-state parameters of hep iron between non spin-polarized calculations and exper- 
iment is significantly larger than for the other two metals studied here (Table ||) or other transition metalsOO. The 
zero pressure volume is underestimated by ~ 9%, and the zero pressure bulk modulus is too stiff by 75% (Table |). 
Especially at low pressure the non-magnetic equation-of-state deviates considerably from experimental values, while 
at high pressure the agreement is very good (Fig. ||) . The stabilization of antiferromagnetic states at low pressure can 
account for some of the discrepancy. For afmll magnetism persists to volumes smaller than Vo, resulting in a larger 
zero pressure volume, reducing the difference with experiment to 5%, and lowering the bulk modulus considerably 
(Table Q). This is still larger than the difference in Vq for cobalt and for cubic iron phases (< 3%)E3'Eil. We attribute 
the remaining discrepancy between low pressure experimental data and the afmll equation-of-state (Fig. |2|) to the 
approximations in GGA and the possible stabilization of more complex spin arrangements than those considered here. 



Elasticity 

The agreement of the calculated elastic constants for cobalt and rhenium with zero pressure experimental results^ 
is excellent with a root mean square error of better than 20 GPa for both metals and both exchange-correlation 
potentials (Fig. [| Tables || andfl]). The initial pressure derivative of the elastic constants for rhenium is also well 
reproduced by the calculations (Fig. |J). LDA and GGA exchange-correlation potentials give almost e qual ly good 
agreement, the minor differences arising primarily from differences in the bulk modulus (Tables |, |l[ and Ell). 

Our elastic constant calculations for rhenium and iron do not agree with the results of lattice strain experiments 
(Fig. ||, Tables [II and IV). For rhenium the overall agreement between these experiments and our elastic constants is 



better than for iron. C\\ and C\i agree well over the pressure range of the experiments, while the other longitudinal 
(C33) and off-diagonal constant (C13) differ significantly (Fig. ||, Table |l|)). The shear elastic modulus (C44) shows 
the largest discrepancy of all elastic constants (factor of 1.5). For iron the results of the lattice strain experiments and 
our calculations are in reasonable agreement for the off-diagonal constants only. The longitudinal moduli we obtain 
at 60 Bohr 3 and 50 Bohr 3 (~ 50 GPa and ~ 200 GPa, respectively) are larger by approximately 50%. This is partly 
related to the overestimated bulk modulus in the calculations. The largest discrepancy, as in the case of rhenium, 
occurs in the shear elastic constants (C44 and Cqq). 

Aggregate properties such as the bulk and shear modulus, and the compressional and shear wave velocity are in 
somewhat better agreement between the theoretical results and the lattice strain experiment for both rhenium and 
iron (Figs. |] and 5|). For rhenium, theory and experiment differ by less than 15% in bulk and shear modulus (Fig. 
|j). For iron the discrepancy is considerable at intermediate pressure but becomes smaller with increasing pressure, 
as already seen for the equation-of-state (Figs. § and ||). At ~ 200 GPa the difference in bulk modulus between GGA 
and experiment is less than 5% and the elastic wave velocities differ by ~ 10%. The shear modulus differs by 25% 
even at high pressure. 

For iron the comparison with previous theoretical results gives a more coherent picture. While the_longitudinal 
elastic constants from our calculations are larger by 10-20% compared to TBB and FP-LMTO resultsEJ (Table |Fv|), 
the elastic anisotropy is similar: the pairs of longitudinal, shear, and off-diagonal elastic moduli display similar values. 
For the TB study this is true over the whole pressure range considered, for the FP-LMTO calculations only at low 
pressure; the ratio of the off-diagonal constants (C12/C13) is strongly pressure dependent in that study, varying from 
0.9 at zero pressure to 0.6 at 400 GPa. 
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IV. DISCUSSION 



We find that the elastic anisotropy (eqs. 14 and |lj) is similar for all three metals studied here. The magnitude of 
the anisotropy is 10±2% for the longitudinal anisotropy and A51, and 30±3% for A52 and is nearly independent of 
pressure (Fig. ||). This is consistent with the experimentally observed behavior of other hep transition metals, all of 
which - except for the filled d-shell metals zinc and cadmium - show anisotropy of similar magnitude (Fig. ^|) . 

Theser«sults can be understood by comparison to a hep crystal interacting with central nearest neighbor forces 
(CNNF)c3. For this model the elastic anisotropy is independent of the interatomic potential to lowest order in P/Cn, 
hence the anisotropy is dependent on the symmetry of the crystal only. Born and HuangO have shown that from this 
CNNF model the elastic constants scale as 32:29:11:8:8 for C 33 :Cn:C 12 :C 13 :C AA , yielding A P = 32/29, A S i = 8/9, 
and A S2 = 45/32 (Fig. |). 

The experimentally determined elastic anisotropics of rhenium and hep iron at high pressure from lattice-strain 
measurements differ substantially from our theoretical predictions, previous theoretical calculations, the behavior of all 
other hep transition metals, and the simple CNNF model (Fig. ^). The shear anisotropy in particular is very different 
in the high pressure experiments as compared with all other relevant results. We suggest that this discrepancy 
may arise from assumptions made inr-fche data analysis. In particular, the assumption that the state of stress on 
all crystallographic planes is identicalta. This condition may not be satisfied in a material undergoing anisotropic 
deformation (e. g. dominated by basal slip), behavior that is observed for many hep transition metals. 

Theory shows much better agreement with lattice-strain experiments in terms of the isotropically averaged moduli. 
Even so, the agreement in the case of rhenium is much better than for iron. In this context it is important to point 
out that part of the discrepancy in the case of iron is due to the mutual inconsistency of the experimentally reported 
elastic constants and isotropic moduli. We have found that the elastic constants reported in Ref. [n] do not yield the 
values of K and \i reported in the same paper. The reason for this discrepancy is unknown. 



V. CONCLUSIONS 



The equations-of-state and the elastic constant tensor at zero pressure and under compression for two ambient 
condition hep transition metals, cobalt and rhenium, and for the high pressure phase of iron, hep, are calculated by 
means of the first principles LAPW method. We find a ferromagnetic ground state for cobalt and an antiferromagnetic 
one for iron, with the antiferromagnetic moment vanishing at 60 Bohr 3 . The equations-of-state for the metals are 
in good agreement with experiment, as arc the clastic constants and pressure derivatives of the elastic constants for 
cobalt and rhenium at ambient pressure. 

Elastic constants for iron under high pressure as inferred from lattice-strain experiments differ significantly from our 
theoretical results. Similarly large discrepancies are also found between theory and high pressure static experiments 
on rhenium. The lattice-strain experiments also lead to large values of the shear anisotropy that differ from that of 
all other open shell hep transition metals. Given the excellent agreement of the theoretical elastic constants for cobalt 
and rhenium with experiment at zero pressure, we suggest that a re-examination of the lattice-strain experiments for 
rhenium and iron is warranted. 
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FIG. 1. Magnetic moment per atom within the muffin-tin sphere for the two antiferromagnetic states of iron considered here 
and the ferromagnetic moment for cobalt as a function of volume. 



FIG. 2. Equations-of-state for the hep metals considered. The upper panel compares the GGA non-magnetic (solid line) 
with the afmll structure (dotted line) for iron; Static experimental data is from Ref. g| (open circles) and Ref. ^ (filled circles). 
The lower two figures show the equations-of-state for ferromagnetic cobalt and non-magnetic rhenium, GGA results are shown 
in solid, LDA in dashed curves. The static experimental data for cobalt are from Ref. 12, for rhenium static (open circles) and 
reduced shock wave data (filled circles) are from Ref. u\ and Ref. nH, respectively. 
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FIG. 3. The elastic constants of hep iron from our calculations are shown in the upper figure. The lines are Eulerian finite 
strain fits to the theoretical results at 45, 50, and 60 Bohr 3 : solid (GGA), dashed (LDA). Lattice strain experiments from Ref. 
and are shown by the open symbols: Cn (o), C33 (A), C12 (v)j Ci3 (^)j an d C44 (□). In the middle panel elastic 
constants for hep cobalt are shown as a function of volume. The curve is again a fit to the calculations at 65 , 70, and 75 Bohr 3 . 
GGA is shown in solid, LDA in dashed lines. At the zero pressure volume they are compared to ultrasonic experiments from 
Ref. |l3| (filled symbols as above). The lower figure shows the equivalent for rhenium with calculations at 85t-93, and 100 Bohr 3 . 
The thick dotted lines indicate the initial pressure derivatives as determined from ultrasonic measurementstS. For lattice strain 
experiments from Ref. [li] open symbols are used again. 

FIG. 4. Isotropic properties for hep rhenium in comparison to experiments. The lower panel shows the bulk (K) and shear 
modulus (fi) of our calculations (GGA) in solid lines. The ultrasonic experiments at ambient condition from Ref. Il5| are shown 
in filled circles with the initial pressure dependence in thick dotted lines. Lattice strain experiments from Ref. Ware shown 
in open symbols. The upper panel uses the same symbols as the lower one for the compressional (v p ) and shear wave velocity 

(«.)• 

FIG. 5. Bulk properties for hep iron in comparison to experiments and the Earth's inner core. The lower panel shows the 
bulk (K) and shear modulus (/x) of our calculations (GGA) in solid lines. Diamond anvil cell experimental results are from 
Refs. |Io| (•) and [l^ (o and □, denoting two different approaches). Ultrasonic measurements in a multianvil experiment (O) 
are from Ref. [jlj as well. The crosses display seismic observations of the inner core. The lower figure uses the same symbols as 
the upper one for the compressional (v p ) and shear wave velocity {v 3 ). 

FIG. 6. As a measure of anisotropy the elastic constant ratios C11/C33, (Cn + C33 — 2Ci3)/4C44, and C44/C66, which govern 
the compressional (Ap) and shear wave anisotropy (Agi and As2, respectively) of the single crystal, are shown as a function 
of the number of d-electrons. The upper figure shows the shear elastic anisotropy Agi, the middle As2, and the lower the ratio 
of the longitudinal elastic constants Ap. For all transition metals crystallizing in the hep phase filled circles are used. The 
dashed lines show the CNNF model predictions. High pressure lattice strain results for iron from Ref. |l(] and for rhenium from 
Ref. ^ are displayed with gray squares. Our results are the open diamonds with the pressure dependence shown in solid lines 
connected to the symbols. 



TABLE I. Equation-of-state parameters from a third order finite Eulerian strain expansion of the energy-volume relation for 
the hep transition metals. Vo, Ko, are the zero pressure volume and bulk modulus, respectively; K the pressure derivative of 
the bulk modulus. For experimental values the bulk modulus is calculated from the elastic constants at ambient pressure. 
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TABLE II. Elastic constants of hep cobalt from theory (GGA, LDA) and experiment. Cee = | (Cn — C12) is added for 
comparison with C44. 
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TABLE III. Elastic constants of hep rhenium from theory (present work: GGA, LDA) and experiment. 6*66 = \ (C11 — C12) 
is added for comparison with C44. 
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TABLE IV. Elastic constants for non-magnetic hep Fe under compression (present work: GGA, LDA); the pressure range 
covered corresponds to approximately 50 GPa to 350 GPa, almost the pressure in the Earth's center. For comparison results 
of other studies at ~ 60 Bohr 3 are included. Cee = \ (Cn — C12) is added for comparison with C44. 
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